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SUMMARY 

Accurate  calculations  to  evaluate  tlxa  performance  of  the 
stoichion-'etric  carbon-oxygen  propellant  syutena  have  been  car- 
ried out  for  nozzle  flow  with  and  without  chemical  reactions  and 
with  and  without  vibrational  adjustment.     The  calculations  show 
that,  for  frozen  chemical  flow,  a  lag  of  vibrational  energy  states 
at  chamber  conditions  nearly  doubles  the  reduction  in  Isp,  as 
comipared  with  flow  in  which  complete  vibrational  equilibrium  is 
maintained.     On  the  other  hand,  lags  in  vibrational  adjustment 
have  practically  no  effect  on  tlic  theoretical  performance  of  hot 
propellant  systems  if  chemical  equilibrium  is  maintained  during 
nozzle  flow.     The  preceding  conclusions  are  in  agreement  with 
the  results  on  other  propellant  systems  obtained  previously  by 
use  of  an  approximate  evaluation  procedure. 
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I.     E^TTRODUCTION 

It  io  wcll-loiown  tliat  the  overall  perfonnance  of  rocket  pro- 
pellant  eysteme  is  dependent  on  the  extent  to  v/liich  physico-chemical 
changes  occur  during  nozzle  flow.   Thus  it  is  customary  to  report 
specific  impulse  data  for  "equilibrium  flow"  (the  chemical  reactions 
are  fast  enough  to  maintain  thermodynamic  equilibrium  at  all  times) 
and  for  "frozen  flow"  (the  chemical  composition  remains  unchanged 
at  the  concentrations  corresponding  to  thermodynamic  equilibrium 
in  the  combustion  chamber).     Practical  experience  has  shown  tliat 
this  method  of  estimating  reasonable  upper  and  lower  bounds  for  per- 
formance evaluation  is  generally  in  accord  with  cxperin^cntal  results.  ♦ 

Implicit  in  the  customary  propellant  evaluation  procedures  is 
the  assumption  that  connplete  thermodynamic  equilibrium  is  main- 
tained during  flow  with  respect  to  excitation  of  internal  degrees  of 
freedom.     It  is  apparent,  however,  that  sizeable  departures  from 
equilibrium  internal  excitation  may  occur,  for  example,  for  the  vi- 
brational degrees  of  freedom,  particularly  for  non-hydrogen  contain- 
ing propellant  systems  in  small  motors  for  v/hich  the  vibrational  re- 

(1  2) 
laxation  times  may  not  be  short  compared  to  the  residence  time.  ^  ' 


♦  The  usual  simplifying  assurrptions  in  propellant  evaluation,  v/hich 
do  not  involve  consideration  of  rate  processes,  are  the  following: 
Thermodynamic  equilibrium  is  reached  in  the  combustion  chamber 
after  adiabatic  reaction;  expansion  of  the  combustion  products  through 
the  Laval  nozzle  is  adiabatic;  the  products  of  combustion  behave  as 
ideal  gases;  the  adiabatic  expansion  may  be  considered  to  involve 
one -dimensional  flow  (parallel  to  the  nozzle  axis)  of  non-viscous  ideal 
gases;  the  velocity  of  the  gases  at  the  nozzle  entrance  position  is  neg- 
ligibly snTiall  compared  to  tliat  at  the  exit  position.     ¥or  further  details 
concerning  customary  evaluation  techniques  see,  for  exaiTiple,   "Quan- 
titative hJvaluation  of  Rocket  Propellants"  by  S.  S.   Penner,  Annerican 
Journal  of  Physics,  Vol.   20,  1952,  pp.   26-31. 
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Approxiinate  eatimates  of  the  effect  of  vibrational  excitation  on 

(3) 
perfori'nance  of  rocketa  have  been  carried  out  previously.  * 

It  is  tlie  purpose  of  the  present  analysis  to  present  the  first 

accurate  evaluation  of  the  probable  conribined  effects  of  chemical 

and  vibrational  lags  during  nozzle  flow.     Performance  calculations 

have  been  carried  out  for  the  following  special  cases: 

(a)  Complete  thermodynamic  equilibrium  is  maintained 
("equilibrium  flow"). 

(b)  Complete  tliermodynamic  equilibrium  is  miaintained 
with  respect  to  internal  energy  states  but  no  chemical 
reactions  occur  during  noszle  flow  ("frozen  flow"). 

(c)  Complete  thermodynam.ic  equilibrium  is  maintained 
except  that  no  adjustments  of  vibrational  degrees  of 
freedom  occur  during  flow  (chemically-cquilibrium 
and  vibrationally-frozen  flow). 

(d)  Complete  thermodynamic  equilibrium  is  nriaintained 
except  that  no  chemical  reactions  and  no  adjustment 
of  the  vibrational  degrees  of  freedom  occurs  during 
flow  (chemically-  and  vibrationally-frozen  flow). 
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U.     BASIC  equations" 

The  requisite  basic  equationo  for  cases  (a)  to  (d)  have  been 

(4) 
given  elsewhere.  *    '    For  the  present  purposes  it  will  be  convenient 

to  express  the  cheinical  composition  by  the  mole  fractions  for  the 

various  chemical  species.     In  terms  of  mole  fractions  it  is  readily 

shown*^  tliat  equations  (7),  (3),  (11)  and  (12)  of  reference  (4)  become, 

respectively: 

(a)    For  equilibrium  flow 


(b)    For  frozen  flow 


n 


Zee  oc      ,    ,       kc 


l^^oko-"oke]  <"'> 


(c)  For  chejiiically-equilibrimrx  and  vibrationally -frozen  flov/ 

i  M^uf  »  txC^  n  Allll  f  (M^  /M   )  #:,  <M^^  -^^        (Ic) 
2      c   e  e  oe      *     c       e    ksl     k«r  ke        kc  -*       * 

(d)  For  chemically-  euid  vibrationally-frozen  flow 

^c^e  •  ^^^e  '  ^' 'oe  +  ^^  'h<c  ^^ke    "  ^^kc  >  <^> 

In  equations  (la)  to  (Id),  u    represents  the  linear  flow  velocity 
at  tlie  nozzle  exit  position,  N,      and  N,      are  the  equilibrium  ^'lolc  frac- 
tions of  species  k  evaluated  at  the  chamber  exit  and  nozzle  exit  posi- 
tions, respectively;  H  ,     and  11  ,      equal  the  total  enthalpy  (sensible 


*For  the  relatively  simple  case  of  flow  with  frozen  vibrational  energy 
statec,  it  is  not  necessary  to  use  the  analysis  given  in  reference  (4) 
since  equations  (la)  to  (Id)  can  be  written  down  immediately. 
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plus  chemical  enthalpy)  per  rnole  of  species  k  if  complete  internal 
equilibrium  exists  at  the  nozzle  entrance  and  exit  positions «  respec- 
tively; the  superscript  vib  signifies  that  only  tlie  vibrational  enthalpy 
is  involved:  the  subscript  o  denotes  complete  energy  equilibriurr^ 
exists  for  all  constituents;  tlie  superscript  o  denotes  complete  chem- 
ical equilibrium  exists  for  all  constituents. 

From  equations  (la)  to  (Id)  it  is  apparent  that  the  evalxiation 
of  u    and  hence  of  Isp  s  u   /f<;  requires  the  calculation  of  equilibrium 

mole  fractions  N,      and  N,      as  well  as  of  total  and  vibrational  en- 
ke  kc 

thalpies.     Since  the  exit  temperature  is  determined  by  the  assump- 
tion that  the  flow  process  is  isentropic,  it  is  obvious  that  the  exit 
temperature  T   ,  and  hence  the  numerical  values  of  N?   ,  hY     ,  and 
H  .    ,  will  depend  upon  the  assumed  type  of  flow  process. 
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iii.    outline:  of  calculations 

For  the  stoichioznetric  solid  carbon-gaacous  oxygen  pro- 
pellant  system   [l  mole  of  C(c)  reacting  with  1  mole  of  0^(g)^     the 
equilibriuni  mole  fractions  at  any  temperature  znay  be  determined 
conveniently  by  writing  tlie  overall  chemical  reaction  in  the  form 

C  +  02"-*(l-x)  CO^  +  X  CO  +  [(x-y)/2]    O^  +  yO  (2) 

Vie  arrive  at  equation  (2)  by  assuming  that  the  first  chem- 
ical reaction  taking  place  in  the  combustion  chamiber  is  C-l-02~*'C02 
which  is  followed  by  a  dissociation  of  part  of  the  carbon  dioxide  to 
carbon  monoxide  and  oxygen  according  to  the  relation  xCO^—*-  xCO 
-f  (x/2)  O.*     This  decomposition  is  then  followed  by  a  dissociation 
of  part  of  the  molecular  oxygen  to  atomic  oxygen  according  to  the 
relation  (y/2)  O, — ^  yO.     In  eqiiation  (2)  x  is  the  number  of  moles 
of  carbon  dioxide  dissociated  and  y/2  is  the  number  of  moles  of 
molecular  oxygen  dissociated. 

From  equation  (2)  we  have 

n^c(2  +  x  +  y)/2  (3) 

where  n.  s  total  number  of  moles  of  reaction  products. 

For  the  chemical  reactions  under  discussion  two  equilibrium 

constants  are  obtained 

1 
••? 

Ki  «  Po  <P02>  ^^^ 

^2'PCO^PO//<PC02>*^/^0  <5) 
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where  K. ,  Kg  and  IC^  art:  in  the  notation  of  data  compiled  from 
reference  (5)*  and  p.    is  the  partial  pressure  of  species  k. 
For  am  ideal  gas 

Pk  ■  "k  Pt^^  <6> 

where  n,    is  the  number  of  moles  of  species  k  and  p.  is  the  total 
pressure. 

By  use  of  equation  (6)  in  equation  (4)  we  obtain 

1  1 

Kj  .  y  [U-y)/2J'7     [2p^/(2  +  x  +  yip  (7) 

and,  similarly,  equation  (5)  becomes 

1  1 


K2«x 


|[(3c-y)/2]  ^  /(l.x)|  [2p^/(2+x+y)]  ^      (3) 


JECquations  (7)  and  (3)  form  a  system  of  two  equations  in  two 

unknowns,  x  and  y.     By  solving  the  equations  for  the  expression 

1  ' 

[p  (x-y)/(2  +  X  +  yfl  ^  *nd  equating  tlie  results  it  is  readily  shov/n 
that 

x^(p^.K2  +  KjK2)  +  x(3K2  -  Y^<.^'li^  «  0.  (9) 

and 

y  «KjX^/[x(Kj-2K2)  +  2K^       .  (10) 

From  equations  (9)  and  (10),  the  quantities  x  and  y,  and  the 
number  of  moles  of  each  of  the  combustion  products,  can  be  deter - 
mined     provided  p.  and  T   ,  the  adiabatic  flame  temperature,  are 
known. 


« Class  notes  for  JP  200,  California  Institute  of  Technology,  1951. 
"^  *  The  effect  of  lag  in  vibrational  energy  states  on  the  equilibrium 
constants  has  been  neglected.     It  will  be  shown  in  Section  IIIC  that 
the  changes  in  equilibrium  constant  will  not  materially  affect  the 
numerical  values  obtained  in  the  present  calculations. 
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In  the  calculations  for  finding  T    the  following  relations  are 


used 


^t''  f^  'Vc  [^^^okc  -  "ok(298. 16°K)1 


and 


av 


O. 


(11) 
(12) 

(13) 


where  Q^     is  the  heat  reieased  by  the  chemical  reaction,    Q    is 
av  c 

the  heat  absorbed  by  the  combustion  products,  and  ^H.,   is  the  stand- 
ard heat  of  formation  of  species  k. 

With  equations  (9)  to  (13)  it  is  found  by  trial  and  error  that 
the  adiabatic  flame  teznperature  is  3750*^K  for  p.  «  20  atmos.      ' 

The  mole  fractions  as  found  by  use  of  equations  (9)  and  (10) 
and  the  average  molecular  weight  M,  as  calculated  frozn  the  relation 


M  a  ZIn,  M 


k*"k 


(14) 


are  given  at  3750  K  and  p.  =  20  atmoe.  and  at  various  teniperatures 
for  pi  s  1  atmofi.  in  Table  I.  Here  N,  and  Kl  arc  the  mole  fraction 
and  molecular  weight,  respectively  of  species  k. 

TABLE  I.     EQUILIBRIUM  MOLE  FRACTIONS  N^  AND  AVER/vGE 
MOLECULAR  WEIGHTS  M  AS  A  FUNCTION  OF  PRES- 
SURE AND  TEMPERATURE  FOR  THE  STOICHIOMET- 
RIC CARBON-OXYGEN  PROPSLLANT  SYSTEM. 


T,°K 

p.,  atmos. 

«CO, 

^^CO 

% 

N^ 

-v.. 

3750 

' 

C.4443 

0.  1787 

0.  087( 

32.  3 

3250 

1 

■ 

0.4631 

J.  1762 

0.  UOf 

31.39 

3000 

1 

:).  365.: 

0.  1587 

0.  047'. 

34.91 

2750 

1 

0.6lv0 

0. 2430 

J.  1150 

0.0  Kv. 

33.09 
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A,     Equilibriuin  Flow 

The  evaluation  of  the  performance  for  equilibrium  flow  re- 
quires tlie  introduction  of  an  cuiditional  assumption,  namelyi  isen- 
tropic  flow  during  the  expansion  procesc  through  the  nozzle.   This 
condition  can  be  expressed  in  terms  of  the  component  entropies  by 
the  relation 

where  5  .      and  S  ,      denote  the  total  molar  enthalpies  at  the  cham- 
ber and  exit  conditions,  respectively  and  R  is  the  molar  gas  con- 
stant. 

The  exit  temperature,   T   ,  cam  now  be  found  by  computing 
tlie  right-hand  side  of  equation  (15)  for  several  temperatures  and 
comparing  the  results  with  the  value  of  the  left-hand  side.  By  using 
the  values  of  M  and  N,   as  given  in  Table  I  and  the  entropies  taDu- 
lated  in  reference  (5),  the  left-hand  side  of  equation  (15)  is  found  to 

be  67.  65  cal.  /°K-mole  for  T     «  3750*^K  and  p^  s  20  atmos.  and  the 

c  c 

right-hand  side  is  foiuid  to  be  66.43  cal.  /  K-tnole  and  68. 16  cal.  /  K- 
mole  at  2900^K  and  3000°K  respectively  for  p.,  «  1  atsnos.     By  linear 
interpolation  it  is  found  that  T^  «  2970°K. 

By  use  of  equation  (la)  it  is  found  that  AH       «  22.  800  kcal/ 
mole  and  17.437  kcal/rnole  at  2900°i;  and  SOOO^^K,  respectively. 
Linear  interpolation  yields  the  result  All       s  19.  04  kcal /mole  at 


the  relation 


T^  s  2970  K  whence  it  follows  that  u     «  2,200  meters /sec.  From 
e  e 

l.p  =\/g  (16) 
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the  specific  impulse  is  seen  to  have  a  value  of  227  seconds. 

D.     Frozen  Flov/ 

For  frozen  flow  equation  (15)  sinaplifies  to 

since  N,      «  N,      and  !  I     «  M   .     3y  proceeding  as  in  Section  UIA 

we  find 

Tc  «  2U0°K; 


from  equation  (lb) 


AH^     a  17. 12  kcal/mole 
oe 


is  obtained; finally  u     s  2105  meters /sec.   and  Isp  =  215  sec. 

C.     Chemically-Ec^uilibrium,  Vibrationally>Frozen  Flow 

For  chemically- equilibrium  but  vibrationaily-frozen  flov/, 
the  isentropic  relation  which  is  used  to  determine  T    becomes 

F,  ^^kc  t^okc  -  P^lnp^-Rlnl^c>  •  ^K /^^e)         • 

(18) 

X    Z:  N°    [5  ,     -Rlnp   .RlnM°     -  {sY^^  -  ^^^^  )1 
,"7,      '^e  L    oke  '^e  ke      *  ke  kc    'J 

For  diatomic  molecules  the  vibrational  contribution  to  the 

(7) 
entropy  is,  in  good  approximation,^ 

-1 


^•^  =  RJ   u  [  cxp(u)  -  l]       -  In  [l  -  e3q>(-u))  I     . 


Here  u  «  hcco*/kT,  h  represents  Planck's  constant,  c  is  the  velocity 
of  light,  k  the  Boltzmann  constant,  and  u>  *  the  v/ave  number  corres- 
ponding to  an  energy  transition  from  the  ground  vibrational  energy 
level  to  the  first  excited  vibrational  state  for  the  non-rotating  diatom- 
ic molecule.     For  the  linear  triatoraic  molecule  CO,,  the  total  vibra- 
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tional  entropy  is  gi  en  by  the  approximate  relation 
^vib 


(7) 


'^  «  R  21 1  u^  [^exp(u, ).l]      -In  ^  l-exptu. )]  j 


(20) 


where  u.  corresponds  to  the  four  normal  vibration  frequencies  of  CO^. 
For  the  present  calculations  the  numerical  values  of  the  cliaracteris- 
tic  temperatures  ©  =  hew*  /k  which  are  listed  in  Table  II  were  used. 


TABLE  U.  CHARACTERISTIC  TEMPERATURES  FOR  O,,  CO,  AND 
CO2 


Molecule 

''?.■ 

CO 

^^'2 

e,V, 

2228 

3067 

955.6,  955.6,   1938.   3493 

The  thermodynaniic  function  S       ,  as  found  by  use  of  equations  (19) 
and  (20),  is  listed  in  liable  III  for  various  ten-iperatures. 

TABLE  in.  S^*^  (cal.  /°K-mole)  AS  A  FUNCTION  OF  TEMPERA- 
TURE FOR  O2,  CO  AND  CO2   / 


3750°K 

3000'^K 

2900°K 

1600^K 

Cora- 
ponent 

li 

-,vib. 

It 

^vib. 

u 

cvib. 
ft 

u 

^vib. 

°2 

2228 

.5941 

3.050 

0.  7427 

2.628 

0.  7683 

2. 5597 

1.393 

1.482 

CO 

3067 

.8178 

2.442 

1.0223 

2.027 

1.0575 

1.9689 

1.917 

0.  972 

1«=02 

19S8 

.5301 

3.271 

0. 6627 

2.  039 

0.6855 

2.7769 

1.243 

1.677 

E^Oz 

955.6 

.2543 

4.709 

0.31ra5 

4.  270 

0.3295 

4. 2025 

0.597 

3.039 

3  CO, 

3493 

.9314 

2.193 

1.1645 

1.774 

1.2045 

1.7334 

2.184 

0.789 

CO2 

Total 

—  -     14.887 

—      13.153 

---     12.920 

---  8.5^44 
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By  use  of  the  appropriate  values    of  N?  and  M  listed  In 
Table  I  and  the  appropriate  values  of  S        listed  in  Table  III  we 
can  now  proceed  to  find  Tc  fro:n  equation  (IS)  by  a  trial  and  err3r 
method  similar  to  tlie  method  described  previously.     With  the  val- 
ues already  obtained  we  find  the  left-hand  side  of  equation  (18)  to 
be  67.  65  cal/'^K-mole  at  T     =  3750°K,  and  p     a  20  atmos.     The 

C  '^C 

right-hand  side  of  equation  (18)  is  found  to  be  ^7.  557  cal/**K-molc 
and  69.  012  cal/°K-niolc  at  2900°K  and  3000^K,  respectively.  By 
linear  interpolation  Hie  result  T     =  2910*'k  is  obtained. 

For  diatomic  molecules  the  vibrational  contribution  to 
the  enthalpy  is,  in  good  approximation,     ' 


,vib 


H^      a  RTu  [exp(u)-l] 


-1 


(21) 


For  the  linear  triatomic  molecule  CO,  the  total  vibrational  en- 


tlialpy  is  given  by  the  approximate  relation: 

2L.   u^   |^e>qp(uj)  -Ij 


H^^^  »  RT 


(22) 


Using  these  relations  the  values  for  H        (in  cal.  /inole)  listed  in 
Table  IV  are  obtained. 

TABLE  IV.  H"""  as  A  FUNCTION  OF  TEMPERATURE  FOR 
O2.  CO  AND  CO2. 


•  '\ 

(c&l/mole)  at 

iVolecul.' 

3750°K 

30G0°i: 

2900°K 

L   ,'..    '        " 

^2 

5455.8 

4G18. 0 

3329.3 

1462. 3 

CO 

4815.0 

3423.  6 

?242.6 

1049. 1 

23,247 

17.424 

16,869 

7116. S 
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As  has  already  been  noted,  tlie  effect  of  lag  in  the  vibrational 

energy  states  on  the  equilibrlun  constants  has  been  neglected.     In 

order  to  see  just  how  nnuch  oi  an  erx'or  this  assumption  introduces 

we  will  now  proceed  with  the  computation  of  K^  for  T  =  3000^K  by 

taking  the  effect  of  lag  in  the  vibrational  energy  states  into  account. 

The  starting  relations  for  the  reaction  CO-  "*"  CO  +  7  ^2  **  ^^^  temp- 

eraturc  T     a  3000°K  are 
e 

.e  «»,,     T.-   c 


AF'=  s  -  RTlnX  ^  (23) 


and 


from  the  given  chemical  re- 
action at  298. 16®K)       (24) 

where  F  is  the  Gibbs  free  energy  as  defined  by  the  relation 

F  «  H  -  TS,  (25) 

and  the  superscript  e  again  denotes  conditions  at  tlie  exit  temper- 
ature.    The  enthalpy,   H,     ,   and  the  entropy,  S,     ,  are  found  from 
the  relations 

respectively.     By  use  of  the  enthalpies  and  entropies  tabulated  in 
reference  (5)  and  of  equations  (25)  to  (27),   togetlier  witii  the  approp- 
riate values  lor  S        and  11        listed  in  Tables  III  and  IV,   the  data 
listed  in  Table  V  are  obtained. 


T 
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In  Table  V,  AF,  represents, the  standard  free  energy  of  formation. 
,  'rom  coiujTm  (5)  of  Table  V  we  have 


o«. 


AF  (resulting  from  tJie  given  chenjical  reaction  at  298. 16  K) 
= -32.  808  -  (-94.  260)  3  +  61.  452 


ami,  using  equation  (24)( 


,-e 


AF     a  +  6.  08  kcal.  /mole. 


Using  tlais  value  for  AF     in  equation  (23)  it  is  readily  shown  that 

.e 


In  K. 


a: 


X'w  J. 


•  3  -  1.  018 


or 


log  K,  =  -  0.442 


at  T  3  3000^K 


On  the  other  hand,  for  complete  energy  equilibrium, 

log  K^  »  -0.4692  at  T  «  3000^K, 
log  I  2  =  -0.  1091  at  T  s  3250**K, 
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ajod,  vising  an  appropriate  linear  interpolation* 

log  K^  —  -  0.44?,  at  1/T  s  .  3314  x  lO"^  or  T  •  ZOn^K, 

From  the  preceding  results  it  is  evident  that  neglecting  the  effect 
of  lag  in  the  vibrational  energy  states  on  the  equilibrium  constant 
K^  at  T  3  3000^iv  has  introduced  an  error  equivalent  to  17°K.     Rep- 
etition of  the  preceding  calculations  shows  tliat  the  error  in  K. 
will  be  even  smaller.     ThercfDre,  we  arrive  at  the  conclusion 
that  the  error  in  the  mole  fractions  at  Te  s  2910'^iC  is  small  enough 
so  that  the  resultant  Te  for  tliie  case  is  within  15°  of  tlie  value  whicii 
would  be  obtained  if  tlie  effect  of  lag  in  the  vibrational  energy  states 
on  the  equilibrium  constants  had  been  taken  into  account. 

Using  the  calculations  given  for  case  (a)  it  is  foxind  that 

AH°^  a  +22.  26  kcal/rr.ole  at  1     a  2910*^K  from  which  it  follows,  by 
oc  e  ' 

using  appropriate  values  of  N,.  and  M  from  Table  I  and  of  11^  ^  from 

Table  IV  in  equation  (Ic)  tliat  ^  M^u^  e  AII^*^  =  +  13.  72  kcal/mole, 

u     s  2,203  meters/sec,  and  I       s  225  seconds. 
e  sp 

^*     Constant-Composition  and  Vibrationally-Frozen  Flow 

i^'or  constant-composition  and  vibr&tionally-frozen  flow  the 
isentropic  relation,  which  is  used  to  deteriviiie  T   ,  becoznes 


n 

k 


^    N°,[  S^^^  -  S^^^  -  (S;^".  S^^^     =  Rln(pc/pe).  (28) 


With  the  values  of  N,  listed  in  Table  I,  the  appropriate  values  of 
S  given  in  Table  IH,  and  the  entropies  tabulated  in  reference  {ij 
it  is  found  from  equation  (28),  by  use  of  a  trial  and  error  rnetlxod. 
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that  T^  e  1600'^K.     Using  the  values  of  N?.  Ubulated  in  Table  I, 

the  appropriate  values  of  H^^    given  in  Table  IV,  and  the  enthalpies 

tabulated  in  reference  (5),  in  eq>-iation  (Id)  it  is  found  that  i  M  u^  = 

c      c  e 

AH^  a  15. 144  kcal/niole,  u^  s  1,981  meters /sec,  and  I       a  202  seconds. 
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IV.  DISCUSSION  OF  RESULTS 

For  convenience  the  results  of  the  present  analysis  are 
•ummarlzcd  in  Table  VI.     Reference  to  Table  VI  shows  that  the 
effect  of  vibrationally-frozen  flow  is  of  no  consequence  when  chem* 
ical  equilibrium  is  maintained.     In  tliis  case  lag  in  vibrational  ad- 
justment lowers  the  perforznance  by  roughly  one  percent. 

TABLE  VI.     SUMMARY  OF  RESULTS 


Type  of  I  lov/ 

,  '  '    ) 

Jn-./s) 

I&p(eec) 

Chemically-cquilibriut.    .zzr. 
cncrgy-cquilibri'      P-^- 

2970 

2220 

227 

Chem.ically-froi:( 
cnergy-equilibri.;      ...> 

?110 

2105 

215 

Chcinically-equilibriu-rc  c^d 
vibrationally-frozen  flow 

2910 

2203 

225 

Chemically-frozen  ^nd 
vibrationally-frozen  flow 

1600 

1981 

202 

On  the  other  hand,  the  effect  of  vibrationally-frozen  flov/ 
for  frozen  chemical  flow  is  such  as  to  decrease  Isp  approximately 
6  /o.     In  this  connection  it  is  of  interest  to  note  that  lags  in  vibra- 
tional adjustznent,  for  a  given  propellant  system,  become  more 
likely  as  the  nozzle  size  is  reduced  because  a  reduction  in  nozzle 
size  means  increased  cooling  rates.     In  particular,  it  has  been 
postulated  that  frozen- vibrational  flow  occurs  in  the  small  nozzles 
used  in  the  pneumatic  method  for  micasuring  temperatures  in 
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combustion  chambers.      ' 

The  calculations  on  the  carbon-oxygen  propellant  system 
have  been  repeated  by  using  an  approximate  procedure/   '  which 
permits  ready  estimates  of  the  effect  of  vibrational  lag  on  per- 
formance for  propellant  systems  for  which  the  usual  theoretical 
performance  data  are  available.    The  results  of  these  calculations 
are  summarized  in  Table  vn,  together  with  data  for  other  pro- 
pellant systems  which  were  obtained  by  utilizing  the  concept  of  a 
constant  average  vibrational  heat  capacity  during  flow.  *   '    Refer- 
ence to  Table  VII  shows  that  lag  in  vibrational  energy  states  pro- 
duces a  similar  effect  on  performance  for  widely  different  propel- 
lant systeiTis.    A  simple  physical  explanation  for  this  observation 
is  obtained  by  noting  that  frozen  vibrational  flow  reduces  the  ef- 
fective heat  capacity  during  expansion  by  nearly  the  same  extent 
for  chemically -frozen  and  for  chcmically-equilibrium  flow.     On 
the  other  hand»  the  effective  heat  capacity  is  much  larger  for  flow 
in  which  chemical  changes  occur  provided  the  propellant  system 
is  sufficiently  hot.     Hence  the  effect  of  vibrational  lags  will  always 
be  more  pronounced  for  flow  without  chemical  change  than  for 
chcmically-equilibrium  flow.   Furthermore,  the  disparity  will  be 
largest  for  tlie  hottest  propellant  systezns.   The  preceding  remarks 
are  obviously  in  accord  with  the  data  given  in  T  ible  \9I  which  show 
a  large  effect  of  vibrational  lag  on  performance  for  cheinicaliy- 
equilibrium  flow  only  for  the  relatively  cool  RFNA-Aniline  propel- 
lant system. 
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'ABLE  VII.     EFFECT  OF  LAG  IN  VIBRATIONAL  ENIIIRGY  STATF?; 
VI  ORMANCE  FOR  VARIOUS  PROPELLANT 


PropcUant  Syatern 

Isp  (energy  eqailibriuii:  iiov/)^- 
Isfj  (vibrationally  frozen  flow)  foi 

ilquiiibriurr!  Flov/ 

^' rosen  !  io 

::_■.>..     ,.:y.c.. 

i. 

A   .       ^  '  .■  ■  J. 

Carbon- ^Oxygen*  * 

1.01 

1.06 

'^  lydrogen-riuorine*  '''■ " 

i.m  ; 

1.  J.I.;. 

:  iydrogen-'-xy^en'-  * <'■ 

- 

RFNA  -  Aniline*** 

1 .  ^  :•} 

1.056 

*  Accurate  calculations. 

*♦  Approximate  calculations. 

***  Results  given  in  reference  (3). 
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